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We calculate the "naive non-abelianization" (NNA) contributions of the higher- twist Feynman 



diagrams to the large-pr inclusive pion production cross section in proton-proton collisions in the 
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case of the running coupling and frozen coupling approaches. We compare the resummed "naive 



non-abelianization" higher-twist cross sections with the ones obtained in the framework of the 



t — " frozen coupling approach and leading-twist cross section. The structure of infrared renormalon 
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singularities of the higher twist subprocess cross section and it's resummed expression are found. 

We discuss the phenomenological consequences of possible higher-twist contributions to the pion 
/\ ' production in proton-proton collisions in within NNA. 
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I. INTRODUCTION 

The hadronic wave functions in terms of quark and gluon degrees of freedom play an 
important role in the quantum chromodynamics predictions for hadronic processes. If the 
hadronic wave functions were accurately known, then we could calculate the hadronic distri- 
bution amplitude and structure functions for exclusive and inclusive processes in quantum 
chromodynamics (QCD). 

The large-order behavior of a perturbative expansion in gauge theories is inevitably dom- 
inated by the factorial growth of renormalon diagrams [1-4]. In the case of QCD, the coeffi- 
cients of perturbative expansions in the QCD coupling a s can increase dramatically even at 
low orders. This fact, together with the apparent freedom in the choice of renormalization 
scheme and renormalization scales, limits the predictive power of perturbative calculations, 
even in applications involving large momentum transfers, where a s is effectively small. In- 
vestigation of the infrared renormalon effects in various inclusive and exclusive processes is 
one of the most important and interesting problem in the perturbative QCD (pQCD). 

The frozen coupling constant approach in Refs. [5,6,7,8,9] was used for calculation of 
integrals, such as 

J 1-x 

It should be noted that, in pQCD calculations, the argument Q 2 of the running coupling 
constant should be taken equal to the square of the momentum transfer of a hard gluon 
in a corresponding Feynman diagram, in both the renormalization and factorization scale. 
But definng in this way, a s (Q 2 ) suffers from infrared singularities. Therefore, in the soft 
regions as x\ —¥ 0, and X2 — > 0, the integral in (1.1) diverge and we need some regular- 
ization methods for a s (Q 2 ) in these regions for their calculation. It is known that infrared 
renormalons are responsible for factorial growth of coefficients in perturbative series for the 
physical quantities. But, these divergent series can be resummed by means of the Borel 
transformation [1] and the principal value prescription [10], and effects of infrared renor- 
malons can be taken into account by a scale-setting procedure a s (Q 2 ) — ¥ a s (exp(f (Q 2 ))Q 2 ) 
at the one-loop order results. 

In this work we apply the running coupling approach [11] in order to compute the effects 
of the infrared renormalons on the meson production in proton-proton collisions in "naive 
non-abelianization" approximation. This approach was also employed previously [12-14] 



to calculate the inclusive meson production in proton-proton and photon-photon collisions. 
The running coupling approach in "naive non-abelianization" approximation pion electro- 
magnetic form factor was computed in [15]. 

A precise measurement of the inclusive charged pion production cross section at s/s = 
62.4 GeV and ^/s = 200 GeV is important for the proton-proton collisions program at the 
Relativistic Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory. 

Therefore, it will be interesting that the calculation of the higher-twist effects on the 
dependence of the pion wave function in pion production at proton-proton collisions by 
the running coupling constant approach. In this respect, the contribution of the higher- 
twist Feynman diagrams to a pion production cross section in proton-proton collisions has 
been computed by using the the infrared renormalon improved distribution amplitude of the 
pion. Also, higher-twist contributions which are calculated by the running coupling constant 
and frozen coupling constant approaches are been estimated and compared to each other. 
Within this context, this paper is organized as follows: In Sec. [TT1 we provide formulas for 
the calculation of the contribution of the higher twist and leading twist diagrams. In Sec. IIHI 
we present formulas and an analysis of the higher-twist effects on the dependence of the pion 
wave function by the running coupling constant approach. In Sec. HVt we give the numerical 
results for the cross section and discuss the dependence of the cross section on the pion wave 
functions. We present our conclusions in Sec. [V] 

II. THE HIGHER TWIST AND LEADING TWIST CONTRIBUTIONS TO IN- 
CLUSIVE REACTIONS 

The higher-twist Feynman diagrams, which describe the subprocess qi + q 2 — > 7r + (7r~) +7 
for the pion production in the proton-proton collision are shown in Fig.l. The amplitude 
for this subprocess can be found by means of the Brodsky-Lepage formula [16]: 

M(s,i) — / dx\ \ dx 2 5(l - xx- X2)§-K(xi,X2,Q 2 )T H (s,i;xi,X2). (2.1) 
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In Eq.(2.1), Tjj is the sum of the graphs contributing to the hard-scattering part of the 
subprocess. In our calculation, we have neglected the pion and the proton masses. The 
Mandelstam invariant variables for subprocesses q\ + q 2 — > 7r + (7r~) + 7 are defined as 

s = (pi+P2)\ i=( Pl - PlT ) 2 , u=( Pl - Pl f. (2.2) 



One of the important moment in our study is the choice of the pion wave functions 
& n (x, Q 2 ) in Eq.(2.1). In Ref.[16] the authors have calculated the contribution of "bubble 
chain" diagrams to the Brodsky-Lepage evolution kernel V[x, y; a(Q 2 )} in the "naive non- 
abelianization" (NNA) approximation and, as a result, have got new, infrared renormalon 
improved distribution amplitude for the meson in the form: 



*«(x, Q 2 ) = Im[x(1 - x)] 1+a J2 a n (Q 2 )A n Cf/ 2+a (2x - l] 



(2.3) 



n=0 



where C n a (2x — 1) are the Gegenbauer polynomials, A n (a s ) are normalization constants, 
a n (Q 2 ) define the evolution of $„■(:£, Q 2 ) with Q 2 and a = a s (Q 2 )(3o/4:7T. In our calculations 
are normalization constants A n (a s ) are given by the expression 

T(3 + 2a) n\ 3 + la + 2n 
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(2.4) 



\/3r(l + a)r(2 + a) (2 + 2a) n 2 + 2a + n ' 

where T(a) is the Euler gamma function, (a) n is the Pochhammer symbol, (a) n = T(a + 
n)/r(o). It should be noted that normalization constants A n (a s ) given by expression (2.4) 
differ from ones of Ref. [17]. If in the Ref. [17] change a by a + 1, then we obtain expres- 
sion (2.4). We have used three different wave functions: the asymptotic (asy), and the 
Gosdzinsky-Kivel wave functions [17]. 



$GKl(x,l4) = ^asy(x,a) 



a 



3/2+a 




(2x - 1) + l.7A 2 C 2 /z+a (2x - 1 



3/2+a, 



®GK2(x,l4) = ^asy(x,a) 



C 3 /2+a (2x - 1) + 1.7A 2 C 2 /2+a (2x - 1) + l.QA 4 C 3 4 /2+a (2x - 1) 



*. i \ r T(4 + 2a) 11+ 

$asy{x, a) = U ", n , \ [x(l ~ X)] 1+a , 



2r 2 (2 + a 
3/2+a ^ , 3 + 2a 



(2.5) 



Cl /2+a (2x - 1) = 1, C' 2 v ^ a (2x - 1) 
,:V2- fVo „ 1N (3 + 2a)(5 + 2a 



((5 + 2a)(2x-l) 2 -l), 



Cf +a (2a; - I) 
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[(9 + 2a)(7 + 2a)(2x - l) 4 - 6(7 + 2a)(2x - l) 2 + 3] 



where /„- is the pion decay constant. The evolution of the wave function on the factorization 
scale Q 2 is governed by the functions a n (Q 2 ), 

"«.(£) lk (x) + C kk (x) 



where 



a n (Q ) = a n (/j, )exp 
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C kk (a s ) = {ot s ^Y{4){- + a + 2k)- ^(- + a)}, 



(2.6) 



(2.7) 
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r . (l + a) 2 r(2a + 4) 



a + l)(a + 2) 2(a + 2) 



'3(2 + a)r(2 + a) 3 r(l -a) { (a + 1 + n)(a + n + 2) 



a + 1 



(a+2+n)-^(a+2)]}. 
(2.8) 



In the limit a — ¥ from infrared renormalon improved distribution amplitude we can obtain 
ordinary distribution amplitude. 

The cross section for the higher-twist subprocess q\q 2 — > 7r + (71-^)7 is given by the expres- 



sion 



where 
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D(t, it) = e\t I dx 
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+ e 2 -u / dx 
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'a a (Ql)$*(x,Ql) 



1 — X 



(2.9) 



(2.10) 



Here Q\ = (x — l)u, and Q\ = —xi, represent the momentum squared carried by the hard 
gluon in Fig.l, e\{e<i) is the charge of qi(q 2 ) and Cp = |. The higher-twist contribution to 
the large-pr pion production cross section in the process pp — > 7t + (tt~) + 7 + X is [18,19]: 
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dx 1 dx 2 G qi / hl (x 1 )G q2 / h2 (x 2 ) — ? ;(qq-nr'y)6(s + i + u). (2.11) 

it at 



In the numerical calculations we denote the higher-twist cross section obtained using the 
frozen coupling constant approximation by (S^ T )°. 

One of the essential problem, extracting the higher-twist corrections to the pion produc- 
tion cross section and a comparison of higher-twist corrections with leading-twist contribu- 
tions. We take two leading-twist subprocesses for the pion production: (1) quark- ant iquark 
annihilation qq — > #7, in which g — > 7t + (tt~) and (2) quark-gluon fusion, qg — > 37, with 
subsequent fragmentation of the final quark into a meson, q — > 7t + (tt~). 



III. HIGHER-TWIST MECHANISM AND THE ROLE OF INFRARED RENOR- 
MALONS 



In this section, we will calculate the integral (2.10) using the running coupling constant 
approach in the naive non-abelianization approximation and also discuss the problem of 
normalization of the higher-twist process cross section in the context of the same approach. 
As is seen from (2.10), in general, one has to take into account not only the dependence of 
<y(Q 2 ) on the scale Q 2 , but also an evolution of $(#, Q 2 ) with Q 2 . Therefore, it is worth 



noting that, the renormalization scale (argument of a s ) should be equal to Q\ — (x — l)u, 
Q\ = —xt, whereas the factorization scale [Q 2 in § M (x,Q 2 )] is taken independent from x, 
we assume Q 2 = p\. It should be noted in Ref.[17] the authors also noted the existence of 
two kinds of power corrections to cross sections; the infrared renormalon ambiguity arising 
from the loop integration and power corrections from regions as, x — > 0, x — >■ 1. The integral 
(2.10) in the framework of the running coupling approach takes the form 

/ (/ i %) = / \^- x • (3- 1 ) 



o 



The a s (Xfi 2 Ro ) has the infrared singularity at x — ¥ 1, for A = 1 — ion- > 0, for A = x and 
so the integral (3.1) diverges. Hence, the integral (3.1) can be found after regularization 
of a s (\ii 2 Ro ) in these end point regions. Such regularization can be fulfilled with the aid of 
the renormalization group equation that allows us to express the running coupling constant 
a s (XQ 2 ) in terms of a s (Q 2 ). The solution of renormalization group equation for the running 
coupling a = a s /n has the form [10] 



a(\) 
Then for a(\Q 2 ), we get 
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(3.2) 



where t = An /a s (Q 2 )(3 = A/a[5 Q . 

Having inserted Eq.(3.3) into Eq.(2.10) we obtain 

Jo (l-x)(ti + lnA) 

. ( h , [' fMJx(l-x)r«ZZo*n(Q 2 )A n C 3 n /2+a (2x-l) 

e 2 ua s (-t)t 2 ^ dx (i- x )( t2 + lnA) (3 ' 4) 

where t\ = 4n/a s (—u)f3 Q and t 2 = 4:7i/a s (—t)(3 . 

The integral (3.4) as we know, still divergent, which can be defined by existing methods. 
Using the running coupling approach it may be found as a perturbative series in a s (Q 2 ) 
with factorially growing coefficient. Making the change variable as z = In A we obtain 

D(t,u) = e 1 ta s (-u)t 1 dx- — ■ — - + e 2 ua s {-t)t 2 I dx— -— — ■ — - (3.5) 

Jo [l-x){t l + z) J (l-x){t 2 + z) 

In order to calculate (3.5) we will apply the integral representation of l/(t + z) [20,21]. 

1 f°° 

= / e~ {t+z >du (3.6) 

(t + z) Jo 



gives 

D{t,u) = eita a (-u)ti / T r h 

Jo Jo I 1 ~~ x ) 

. ( fu f l r ® M (x,p 2 T )e-^ u dudx 

e 2 ua s {-t)t 2 / / j- r • (3.7) 

Jo Jo I 1 ~ x ) 

In the case <& asy (x,a) for the D(t,u), it is written as 

.* s 47rLeii r(4 + 2a) Z" 00 , , „ . 



4nUe 2 u T(4 + 2a) f°° _ <2UD . o 



(3.* 



and for $G^ 1 (a;, Q 2 ) wave function 

£>(£,«) = lElpl ^ + 2a) ■ r^1Ai(«.)fl(2+a ) l + a-u) + 1.7A 2 («,)^. 
Po 2V3r^(2 + a) Jo ^ 

Att p(^) 7fc (x) + C fcfc (^ V 
exp — - / ax ((5 + 2a)(45(4 + a, 1 + a — u) — 

\0oJ «.(£) X 2 J 

45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - «)) - 5(2 + a, 1 + a - u)] + 

AirLe 2 u T(4 + 2a) [°° , .„. , . ,„. , „ „ , ,3 + 2a 

Po 2V3r z! (2 + a) Jo ^ 

Ayr p(4) 7fc(a;) + Cfcfc(a;) \ 

ea:p — - / ax ((5 + 2a)(45(4 + a,l + a — u) — 

\PoJ a .o%) x 2 J 

45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - u)) - 5(2 + a, 1 + a - u)]. (3.9) 

also for $gx2(^, Q 2 ) wave function 

£>(f,fi) = ^^- ^ 2 " } , • rdue-^[A (a s )B(2 + a,l + a-u) + 1.7A 2 (a s ) 3 ^. 

Po 2V3r^(2 + a) Jo 2. 

exp \^l r {PT) ik(x)+c kk (x) \ + 2a)(4jB(4 + a> x + a _ u) _ 

\PoJ a,Q%) x J 

45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - u)) - 5(2 + a, 1 + a - u)] + 
exp I ^ r J, ' r) 7fcW + C fcfc (x) ^\ + ^ + 2a)(16jB(6 + Q x + a _ u) _ 

\PoJ a.(j%) X 2 J 

325(5 + a, 1 + a - u) + 245(4 + a, 1 + a - u) - 85(3 + a,l + a-u) +5(2 + a, 1 + a - «))- 



6(7 + 2a)(45(4 + a, 1 + a - u)- 
45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - u)) + 35(2 + a, 1 + a - «)]- 

dMe" t2 "[A (a s ) J B(2 + a, 1 + a - u) + 1.7A 2 {a s 



AirUe 2 u T(4 + 2a) [<*> J _.__ t2ttr , ,_. ND/n ,....,, . .. N ,■,,,, ,_ ,3 + 2a 
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4vr p^fcM + CW*) 



exp[— ' v y w cfa| (5 + 2a)(4B(4 + a,l + a-tz)- 

45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - u)) - 5(2 + a,l + a- u)] + 

(3 + 2a)(5 + 2a) 
l.bAAa s ) — 
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325(5 + a, 1 + a - it) + 245(4 + a, 1 + a - it) - 85(3 + a,l + a-u) +5(2 + a, 1 + a - u))- 



Att /• Qs( ^) 7fc ( x ) + C fcfe (x) , \ .. w 

expl — - v ; - v ; cfo ((9 + 2a)(7 + 2a)(165(6 + a, 1 + a - «)- 



6(7 + 2a)(45(4 + a, 1 + a - u)- 

45(3 + a, 1 + a - u) + 5(2 + a, 1 + a - u)) + 35(2 + a, 1 + a - «)]. (3.10) 

where B(a,j3) is Beta function. 

The inverse Borel transformations have the infinite number of infrared renormalon poles 
at the points uq = k + a in the Borel plane. According to Ref. [4] that infrared renormalon 
pole at u = k correspond to a power-suppressed correction as (A 2 /Q 2 ) k to a physical 
quantity under consideration. If the renormalon pole is located ai Uq = k + a then its 
contribution is order (A 2 /Q 2 ) k /e or (A 2 /Q 2 ) k /e 2 . Therefore, our expressions (3.8)-(3.10) 
takes into account the power-suppressed corrections Ck(Q 2 )(A 2 /Q 2 ) k , k = 1,2,3,.. to the 
inclusive pion production cross section in proton-proton collisions. The coefficients Ck{Q 2 ) 
of these corrections also depend on the chosen pion wave functions. It should be noted 
that here we neglect infrared ambiguities SCk(Q 2 ){A 2 / Q 2 ) k producing by the principal value 
prescription itself, which have to be canceled by ultraviolet renormalon ambiguities of higher 
twist corrections to cross section and do not estimate 5Ck{Q 2 )- The structure of the infrared 
renormalon poles in Eqs. (3.8-3. 10) strongly depends on the wave functions of the pion. In 
the numerical calculations we denote resummed higher-twist cross section by (S^ T ) res . 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section, the numerical results for the "naive non-abelianization" contribution 
of higher-twist effects to large-p^ inclusive pion production cross section in the process 
pp — > 7r + (or 7r~)7 + X in the case higher-twist contributions calculated in the context of 
the running coupling and frozen coupling approaches on the dependence of the infrared 
renormalon improved distribution amplitude of the pion are discussed. In the numerical 
calculations for the quark distribution function inside the proton, the MSTW distribution 
function has been used [22]. The gluon and quark fragmentation functions into a pion 
has been taken from [23]. The results of our numerical calculations are plotted in Fig.2- 
Fig.9. In Fig. 2 - Fig. 4 we show the dependence of the higher-twist cross sections (S^ 1 ) , 
(£^ r ) res calculated in the context of the frozen and running couplings constant approaches 
and the ratios R = (Sf + T ) re VSf + T )°, (Sf + T )°/(S^) (Sf + T ) re V(S^) as a function of the 
pion transverse momentum pt for different pion wave functions at y — 0. It is seen from 
Fig. 2 that the higher-twist cross section is monotonically decreasing with an increase in the 
transverse momentum of the pion. In Fig.3 - Fig. 4, we show the dependence of the ratios 
R = (£f + T ) re 7(£f + T )°, (Sf + T )°/S^, (£f + T ) re 7£^ as a function of the pion transverse 
momentum px for different pion wave functions. Here S^+ is the leading-twist cross section, 
respectively. As shown in Fig.3, in the region 2 GeV/c < pr < 4 GeV/c resummed higher- 
twist cross section is suppress by about 1-2 orders of magnitude relative to the higher-twist 
cross section calculated in the framework of the frozen coupling approach. In Fig. 4, we 
show the dependence of the ratios (S^j r )°/E^+, and (S^J r ) res /S^+ as a function of the 
pion transverse momentum px for different pion wave functions. It is observed from Fig.4 
that, the ratio (E^T)°/E^+ with increasing transverse momentum of pion decrease and 
has a minimum approximately at the point pt = lOGeV/c. After that, the ratio increase 
with increasing in the pt transverse momentum of the pion. Also, as shown in Fig.4, the 
ratio (S^ r ) res /S^+ decrease with increasing in the pt transverse momentum of the pion. 
In Fig.5, we have depicted higher-twist cross sections (X^) ; (Y,^+) res , as a function of 
the rapidity y of the pion at y/s = 62.4 GeV and px = 4.9 GeV/c. Figure show that 
higher-twist cross section in case frozen and running coupling constant approaches have 
a different distinctive. In the region (—2.52 < y < 1.22) the higher-twist cross section 
{T l ^+)° for all wave functions increase with an increase of the y rapidity of the pion and 
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has a maximum approximately at the point y = 1.22. But, the resummed higher-twist 
cross section (S^ r ) res for all wave functions increase with an increase of the y rapidity 
of the pion and has a maximum approximately at the point y — — 1.92. As is seen from 
Fig. 5 in the region —2.52 < y < 2 resummed higher-twist cross section for <&gki(x, Q 2 ) 
and §gk2{ x , Q 2 ) is suppressed by about one order of magnitude relative to the resummed 
higher-twist cross section for Q asy (x,a). As is seen from Fig.6 ratio R = (£^j r ) res /(£^4 r )°, 
for all wave functions increase with an increase of the y rapidity of the pion and has a 
maximum approximately at the point y = —1.92. In the region —2.52 < y < —1.92 
resummed higher-twist cross section is suppress by about two orders of magnitude relative 
to the higher-twist cross section calculated in the framework of the frozen coupling approach. 
Besides that, the ratio decreases with an increase in the y rapidity of the pion. Analysis 
of our calculations shows that (S^ 1 ) , (Y,^+) res higher-twist cross sections are sensitive to 
the choice of the infrared renormalon improved distribution amplitude of the pion. We have 
also carried out comparative calculations in the center-of-mass energy y'i = 200 GeV and 
obtained results are displayed in Fig.7-Fig9. Analysis of our calculations at the center-of- 
mass energies \fs = 62.4 GeV and ^/s = 200 GeV, show that the increasing in the beam 
energy contributions of higher twist effects to the cross section decrease by about one-two 
order. In our calculations of the higher-twist cross section of the process the dependence 
of the transverse momentum of pion appears in the range of (10 -8 — 10~ w )mb/ GeV 2 . We 
think, that higher-twist cross section obtained in this work should be observable at RHIC. 

V. CONCLUSIONS 

Proton-proton collisions are known to be the most elementary interactions and form the 
very basis of our knowledge about the nature of high energy collisions in general. Physicists, 
by and large, hold the view quite firmly that the perturbative quantum-chromodynamics 
provides a general framework for the studies on high energy particle-particle collisions. Obvi- 
ously, the unprecedented high energies attained at Large Hadron Collider offer new window 
and opportunities to test the proposed QCD dynamics with its pros and cons. However, we 
should remember that LHC opens a new kinematical regime at high energy, where several 
questions related to the description of the high-energy regime of the QCD. Consequently, 
studies of proton-proton interactions at the RHIC and LHC could provide valuable infor- 
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mation on the QCD dynamics at high energies. In this work the " naive non-abelianization" 
(NNA) contributions of the higher-twist Feynman diagrams to the large-pr inclusive pion 
production cross section in proton-proton collisions are calculated. For calculation of the 
higher-twist cross section the running coupling constant approach is applied and infrared 
renormalon poles in the cross section expression are revealed. Infrared renormalon induced 
divergences is regularized by the means of the principal value prescripton and the resummed 
expression (the Borel sum) for the higher-twist cross section is find. It is abserved that, 
the resummed higher-twist cross section differs from that found using the frozen coupling 
approach, in some regions, considerably. Also we have demonstrated that higher-twist con- 
tributions to pion production cross section in the proton-proton collisions have important 
phenomenological consequences. Future RHIC and LHC measurements will provide further 
tests of the dynamics of large pt hadron production beyond leading twist. 
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FIG. 1: Feynman diagrams for the higher-twist subprocess, q\q2 —> vr + (or tt ) / y. 
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FIG. 2: Higher-twist 7r + production cross section (S + ) as a function of the pr transverse mo- 
mentum of the pion at the cm. energy ^/s = 62.4 GeV. 
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FIG. 3: Ratio R = (S + ) res /{Y, + )°, where higher-twist contribution are calculated for the pion 
rapidity y = at the cm. energy -y/i = 62.4 GeV as a function of the pion transverse momentum, 
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FIG. 4: Ratio (Yi H +)/Y, L +, as a function of the px transverse momentum of the pion at the 
cm. energy y/s = 62.4 GeV. 
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FIG. 5: Higher-twist 7r + production cross section (E + ) , as a function of the y rapidity of the 
pion at the transverse momentum of the pion pt = 4.9 GeV/c, at the cm. energy \fs = 62.4 GeV . 
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FIG. 6: Ratio R = (T.^+) res /(T l ^) , as a function of the y rapidity of the pion at the transverse 
momentum of the pion px = 4.9 GeV/c, at the cm. energy ^fs = 62.4 GeV . 
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FIG. 7: Higher-twist n + production cross section (Y, H +) as a function of the pt transverse mo- 
mentum of the pion at the cm. energy y/s = 200 GeV. 



16 



10" 1 



10"' 



o 
ii 

— 10" 3 



w 
sp 10" 4 

10" 5 



10"' 



\ 



20 



■ asy(0) 
GK1(0) 
GK2(0) 
asy(res) 
GK1(res) 
GK2(res) 



40 60 

p T , GeV/c 



80 



100 



FIG. 8: Ratio (E + )/E +, as a function of the pr transverse momentum of the pion at the 
c.m.energy y^ = 200 GeV. 




FIG. 9: Ratio R = (E + ) res /(E + )°, as a function of the y rapidity of the pion at the transverse 
momentum of the pion pr = 15.5 GeV/c, at the cm. energy y/s = 200 GeV. 
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